Propagation and Backscattering of Mechanical Impulses in a Gravitationally Loaded 
Chain: Dynamical Studies and Toy Model Based Phenomenology 



in 
o 
o 

(N 
-i— > 

o 
O 
o 

(N 



c3 



i 

c 

o 
o 



in 
in 
o 

in 
o 

c3 



i 

C 

o 
o 



X 



Edgar AvalosQ T. R. Krishna MohanQ and Surajit Ser|3 

Department of Physics, State University of New York. Buffalo, New York 1^260-1500 

(Dated: February 2, 2008) 

We recently introduced a simple toy model to describe energy propagation and backscattering in 
complex layered media (T.R. Krishna Mohan and S. Sen, Phys. Rev. E 67, 060301(R) (2003)). The 
model provides good qualitative description of energy propagation and backscattering in real soils. 
Here we present a dynamical study of energy propagation and backscattering in a gravitationally 
loaded granular chain and compare our results with those obtained using the toy model. The 
propagation is ballistic for low g values and acquires characteristics of acoustic propagation as g is 
increased. We focus on the dynamics of the surface grain and examine the backscattered energy 
at the surface. As we shall see, excellent agreement between the two models is achieved when we 
consider the simultaneous presence of acoustic and nonlinear behavior in the toy model. Our study 
serves as a first step towards using the toy model to describe impulse propagation in gravitationally 
loaded soils. 

PACS numbers: 46.40.Cd,45.70.-n,43.25.+y 



I. INTRODUCTION 

Acoustic imaging of buried objects in a non-linear 
medium like nominally dry soil is still an open prob- 
lem 0. Amongst its main applications, we can mention 
the problem of locating antipersonnel land mines, human 
remains for forensic investigations, hidden underground 
structures of archaeological importance etc. It has been 
shown that gentle mechanical impulses 2] can be used to 
detect buried objects at depths of a meter or so in nomi- 
nally dry sand beds. Imaging requires the understanding 
of pulse propagation in nonlinear media. Besides the un- 
derlying non-linearity due to the Hertzian contact forces 
between the grains (see, for example, Q), gravitational 
loading is an added feature that affects the dynamics. 
Before we can analyze the gravitational effects in 3D sys- 
tems, it is reasonable first to look at the problem in ID. 
Further motivation for ID studies can be found in 

Impulse propagation in a gravitationally loaded chain 
with N grains and a perfectly reflecting boundary at the 
bottom has been studied lEJaQl. When 5 = 0, one finds 
solitary wave propagation [5|, Q as shown via velocity vs. 
position plots made at different times in Figs, ^a)-(c) 
(see later for simulation details). In all the panels of 
Fig. 2] the numbers in the bottom left corner indicate 
the times at which the snapshots have been taken and, 
in brackets, the round trip time of the pulse is given. 

The velocity of the solitary wave is a function of its 
amplitude 0, H, 0, |HJ • Each solitary wave is reflected 
at the bottom boundary. During reflection, the solitary 
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wave forms secondary solitary waves ^2 Ell ■ Fur- 
ther, secondary solitary waves are generated from col- 
lisions with other solitary waves Q. The formation of 
secondary solitary waves progressively reduces the am- 
plitude of the solitary waves in the system. In Fig.^c), 
we find that the amplitude of the primary solitary wave 
has diminished considerably after 2.6 roundtrips. 

When g > 0, one finds a solitary wave like front with a 
tail that elongates with time E 13. In Figs.UId)- 
(/), we have shown snapshots from the propagation of the 
impulse in a system with g — 0.1 and with open bound- 
ary condition at the top. We have the same boundary 
conditions in Figs. PJg)-(t) with g — 1.0. We see that 
the impulse propagates faster with increased g. 

The scheme of this paper is as follows. Section II sur- 
veys briefly the 1-D model with the simulation details. 
Our focus is to probe impulse propagation and backscat- 
tering by simply considering the dynamics of the surface 
grain. In Section III, we construct a version of the toy 
model and recover the results obtained in section II. Sec- 
tion IV summarizes the conclusions for the ID study and 
assesses the usefulness of this study for further 3D anal- 
ysis. 



II. SIMULATION MODEL 

We model the granular chain as a collection of 10000 
spherical grains which are placed in contact with one an- 
other and loaded by a gravitational field. The interaction 
between every pair of spheres of radius R is driven by 
Hertz law 013, i-e-, it is assumed that spheres labeled 
as i and i + 1 are interacting with a potential proportional 
to the overlaps, in the contact region, 
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FIG. 1: In Fig. n a )~( c )i .9 = and we see solitary wave 
propagation. In Fig. ffl d)-(f). g = 0.1 and in Fig. nTq'l- 
(i), g = 1.0. In both cases, open boundary conditions are 
employed at the top with closed boundary conditions at the 
bottom. We see solitary wave like front with a long tail for 
non-zero g. The text in the lower left hand corner indicates 
the times at which the snapshots are taken and, in brackets, 
is given the roundtrip time, i.e. the time taken by the pulse 
to travel up and down the chain once. The y-axis ranges have 
been marked in the leftmost panel for each row of panels while 
the x-axis range is the same in all and labelled in the top left 
panel, Fig.^a), where j is the position of the particle in the 
chain and Vj its velocity. 



where a is a constant that depends upon the material 
characteristics of the grains, r^+i is the separation dis- 
tance between the centers of the grains i and i + 1 and 
<5^i_|_i = 2R — rij + i is the grain overlap. The equation of 
motion for grain i can be written as [iff, 
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where m and g are the mass of the grain and the value of 
gravity, respectively. The system dynamics is obtained 
by time integration of the coupled Newtonian equations 
of motion via the velocity- Verlet algorithm [ijj]. We ig- 
nore the role of dissipation in the present study; dissipa- 
tion effects can be built in later [18j . We set 2R and m to 
1, and a is set equal to 5657. We find that an integration 
time step of 1 x 10 -5 provides excellent enrgy conserva- 
tion; decreasing the time step further does not improve 




FIG. 2: In Fig. g = 0, in Fig. g = 0.1, and in 

Figs. I3c)-(d), g = 1.0. Figs. I3a)-(c) have closed boundary 
conditions at the top while Fig. [21 d) has open boundary con- 
dition at the top. Typical integrated backscattered energy, 
E bs (t), plots are shown in Figs.|5Je)-(/) where Fig.^e) is for 
g = 0.0 and Fig. |2I/) for g = 0.1 with top closed boundary 
condition (cf. Fig. The center of mass, fi, oscillations 

typically found in these systems are given in the third row 
of panels. Fig. 121c?) is for g — 0.0 while Fig. \Rh) is for the 
top closed boundary condition (g = 0.1; cf. Fig. and 
Fig- EI/))- Fig- 1210 i s f° r t ne top open boundary condition 
(cf. Fig. 121^))- The x-axis remains the same in all panels 
as labelled in Fig. |2[a) with the time, t, given in arbitrary 
units (a.u.) of the simulation. j/-axis ranges for the first row 
is given in the leftmost panel, and separately labelled for each 
panel in the bottom row. For the middle row, the first panel 
on the left has a different y-axis range as shown and, for the 
last two panels, it is as shown for the middle panel. 



the accuracy of our results. Particles are labeled start- 
ing from the bottom so that the iV tn particle is at the 
top. While the bottom boundary condition has been kept 
perfectly reflecting in all the simulations, both open and 
closed boundary conditions have been employed at the 
top in different simulations to investigate the differences 
induced by the particular choice. 

The dynamics is initiated via an initial velocity per- 
turbation at the top of the chain, vif(t = 0) = Vi n , 
with viit = 0) = for < i < N; we have employed 
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Vin = —0.01 in the studies reported here. We subse- 
quently monitor the backscattered energy received at the 
surface, e bs (t) = v%(t)/2. The time integrated backscat- 
tered energy, E bs (t) , has also been computed by summing 
the entire sequence of energy packets received at the sur- 
face, i.e. E bs (t) = £t>=o e bs (t'). 

In Fig.Eto-d), we show plots of e bs (t) vs. t. Fig.H» 
is obtained with g = 0.0. We have two prominent peaks 
at early times followed by a continuous distribution of 
much smaller peaks. Clearly, the smaller peaks originate 
in the breakdown of the solitary waves at the boundaries 
and through mutual interactions (this issue is discussed 
in detail in 01). The two large peaks at early times are 
separated by roundtrip times of the system. Fig. |2[&) 
corresponds to g = 0.1, with reflecting boundary con- 
ditions at the top in which the surface grain is only al- 
lowed to move in one direction into the chain. In this 
case, there are large peaks separated by roundtrip times 
of the system, but there is much activity in between them 
indicating the presence of acoustic like rattling through- 
out the system. Increasing g to 1.0 (Fig. Etc)), with the 
same boundary conditions, does not reveal any noticeable 
change in the pattern except for increasing the density 
of peaks, caused by the smaller roundtrip times. If wc 
keep g = 1.0 and change to open boundary conditions, 
meaning that the surface grain is allowed to move up and 
down, we get the pattern shown in Fig. ^d) where we 
notice that the number of larger peaks are reduced. We 
attribute this to increased energy sharing caused by the 
open boundary at the top. 

The E bs (t) vs. t plots show a distinct difference be- 
tween the ballistic and acoustic type propagation cases. 
While Fig. |^e) is for the g = case, Fig. EI/) i s typical 
of g > cases (here, g = 0.1 with closed boundary at 
the top). The increase is linear for non-zero g whereas 
it only approaches linear for g = 0.0 in the later stages. 
The steps in E bs (t) indicate the arrival of peaks with 
the quiescent period in between marked by plateaus; the 
length of the plateau indicates the time period between 
peaks. 

The pattern of arrival of large peaks introduce modula- 
tions in the amplitudes of e bs {t). Nevertheless, there are 
no simple patterns in the distribution of peaks of e bs (t) 
vs. t in Figs. [21 a-d). We note that the boundary con- 
ditions affect e bs (t). A possible modulating mechanism 
could be due to the center of mass (com) oscillations of 
the system. We show these oscillations (the "breathing" 
of the chain) in the third row of Fig. [21 for three typical 
cases. With 5 = (Fig. Efg)), the significant com oscil- 
lations appear only after the solitary waves have broken 
down. For non-zero g, the com oscillations are affected 
by the surface boundary. Fig. [!J h) is for g = 0.1 and with 
closed boundary conditions at the top while Fig. [21 i) is 
typical of open boundary conditions at the top (in this 
case, g = 1.0). We now turn to the toy model to see 
whether these patterns in e bs (t) can be reproduced. 



III. COMPARISON WITH TOY MODEL 



We had introduced a toy model in an earlier paper 
where the energy propagation in a vertical alignment of 
masses was considered in a simplified manner |4|. At 
time t — 0, we set initial energy E = 1 for layer one 
and zero for the rest. At t = 1, the first layer in the 
vertical chain transfers p (< 1) of the impulse energy to 
the second layer, and retains (1— p). At subsequent times, 
the impulse will propagate in the same fashion all the 
way down the chain. After each transfer of energy, the 
phase of the mass reverses so that it will interact with its 
adjacent layer in the opposite direction in the following 
time step. The interaction between two adjacent layers 
occurs in the following two ways: (i) equipartition case, 
and (ii) exchange case. In the equipartition case, the two 
interacting layers will come away from the interaction 
with equal amounts of energy; we add up the individual 
energies of the two layers and divide the sum equally 
between them. In the exchange case, we let the layers 
exchange their energies; the two interacting layers, after 
the interaction, come away with the energy of the other. 

In Fig. [21 we show the results where the model has, for 
the first time, both equipartition and exchange. Our goal 
is to develop the propagation patterns seen in Figs. [21 b- 
d). The model is based on the exchange case, but we 
allow equipartition centers to develop, as time goes by, at 
different points in the chain. In Fig. [21 a), we have shown 
a case where a single equipartition center is kept fixed at 
layer number 200 in a chain of 500 layers. Fig. [21 6) is for 
a similar case but with the equipartition center located at 
layer number 350. We see that the patterns in e bs (t) do 
change slightly depending on the location of the isolated 
equipartition center, but, in both cases, it is seen that 
the patterns are broadly similar to the non-zero g cases 
seen in Figs. [21 b)- (d). 

This has prompted us to try averaging over the pat- 
terns resulting from differing locations of the isolated 
equipartition centers and, indeed, the plot of Fig. [21 c) 
shows that such averaging does retain the similarity with 
plots of Figs. Wib)-(d). To obtain this fi gure, we aver- 
aged over, along with the cases given in Fig. [21 a) and 
Fig. [21 b), two other cases with isolated equipartition cen- 
ters, in each case, located at layer numbers 300 and 390 
respectively. If we vary the exchange coefficient p, we 
get similar patterns (not shown here) but with the dif- 
ference that the growth in E bs (t) is proportional to the 
p value. This indicates that p could also be used, along 
with the equipartition centers, to capture the transition 
from ballistic propagation to acoustic propagation. In 
Fig. [21 d), we show results with multiple (three, in this 
cases) equipartition centers in the system, at layer num- 
bers 125, 230 and 390; p = 0.8 for this case. We see 
that E bs (t) grows by as much as in Figs, ^a)-(c) but 
the attenuation in peaks of e bs (t) is very pronounced and 
similar to the pattern seen in Fig. \^d). 
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IV. SUMMARY AND CONCLUSIONS 




FIG. 3: All the panels along with the insets have the axes 
ranges as labelled in the top left panel. Fig. 1ST a) is for the 
case of a single equipartition center in the toy model located 
at particle no. 200 out of a total of 500 particles in the chain, 
and Fig- EI ^) f° r the case, similarly, with equipartition center 
at particle no. 350; the exchange coefficient p = 0.8. Data 
for these two cases have been averaged, along with data for 
cases with single equipartition centers located at 300 and 390 
respectively, and presented in Fig. EI c). In Fig. EI d), we have 
multiple equipartition centers located at particle nos. 125, 
230 and 390 with p = 0.8. 



In this paper, we have studied the problem of impulse 
propagation in a gravitationally loaded granular chain. 
We have focused our attention on the backscattered en- 
ergy received at the surface after an impulse has been 
initiated into the system. Our studies have been car- 
ried out using two different approaches — (i) using New- 
tonian dynamics in a non-disspative system to describe 
backscattering at the surface of the system and, (ii) using 
the toy model ]ij to recover the behavior in (i). 

Our results show that the toy model is capable of re- 
producing the correct form of the backscattered energy 
as a function of time in a gravitationally loaded chain. 
Gravitational loading introduces acoustic-like oscillations 
in the system. Such oscillations vanish and the system 
ends up propagating solitary waves when gravitational 
loading is zero. To achieve these descriptions, one must 
generalize the earlier version of the toy model and 
incorporate both exchange (represents non-linear) and 
equipartition (represents acoustic) effects. We find that 
the role of the equipartition (acoustic) effect in the toy 
model is rather dominant and only limited equipartition- 
ing of energy gives the correct backscattering behavior. 
Our studies confirm that mechanical propagation in a 
granular chain is strongly nonlinear even in the presence 
of gravitational loading. 
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